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ABSTRACT 


Lambert  developed  equations  relating  times  of  transit 
between  two  points  in  space  and  the  semimajor  axis  of  conics 
passing  through  these  two  points  when  the  two  radii  and  the 
chord  are  given.  Special  types  of  problems  can  often  best  be 
solved  by  alternate  methods  that  have  been  developed,  but 
for  a  general  study  of  connecting  two  points  in  space  with  a 
conic  section,  with  no  special  constraints  other  than  time, 
Lambert's  equations  seem  to  be  best  suited.  This  paper 
represents  an  expository  summary  of  the  mathematical 
methods  and  techniques  involved. 
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I.  INTRODUCTION 

Statement  of  Problem 

Given  two  positions  in  space  and  their  distances  from 
a  central  body,  find  an  equation  that  relates  time  of  transit 
with  a  geometrical  orbital  parameter;  namely  the  semi¬ 
major  axis.  The  solution  to  this  problem  is  of  interest  in 
orbit  determination  since  knowledge  of  positions  and  times 
are  usually  involved. 

Kepler  derived  an  equation  which  relates  time  from  perifocal  passage 
with  semimajor  axis  and  eccentricity,  given  the  eccentric  anomaly  of  a  point. 
This  equation  is  most  often  used  to  find  the  eccentric  anomaly  when  time, 
semimajor  axis,  and  eccentricity  are  known  (as  is  the  case  in  ephemeris 
prediction).  It  could  be  used  in  a  "direct"  solution,  that  is,  finding  the  time, 
given  the  other  quantities  (as  is  the  case  in  finding  the  time  of  change  of 
phase  in  the  patched-conic  method  of  trajectory  analysis). 

Lambert  studied  the  problem  stated  and  found  that  by  proper  substitu¬ 
tions  in  Kepler's  equations  he  could  eliminate  the  eccentricity  and  so  de¬ 
rive  an  equation  relating  time  and  semimajor  axis  when  two  radii  and  the 
chord  are  given.  Other  methods  of  orbit  determination  have  also  been 
developed,  such  as  the  Lagrange,  Gauss,  and  Gibbs  methods  or  combina¬ 
tions  and  modifications  of  them.  Which  method  to  use  for  a  given  set  of 
known  data  is  a  long  study  in  itself  and  beyond  the  scope  of  this  paper. 
Lambert's  method  is  completely  general  and  the  mathematical  method 
remains  valid  for  problems  involving  long  arcs.  Thus,  for  problems  where 
the  length  of  the  arc  is  unknown  (and  could  be  very  long)  this  method  is 
extremely  useful.  However,  for  certain  types  of  problems  there  is  no  doubt 
that  other  methods  can  often  be  more  efficient. 

II.  BASIC  GEOMETRICAL  RELATIONS 

Given  the  attracting  focus,  there  exist  an  infinite  number  of  conic 
sections  passing  through  two  given  points  in  space.  If  in  addition,  a 
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semimajor  axis  is  given,  Table  1  shows  the  necessary  and  sufficient  condi¬ 
tions  to  classify  the  type  of  resulting  conic  sections. 

Table  1.  Classification  of  Conics. 

If:  Then: 


4a  >  r  ,  +  r  +  c 
d  a 

2  ellipses 

4a  =  r  ,  +  r  +  c 
d  a 

1  ellipse 

4a  <  r  ,  t  r  +  c 
d  a 

no  ellipses 

4a  >  c  -  r  ,  -  r 
d  a 

2  hyperbolas 

(concave  branch) 

4a  <  r  .  +  r  -  c 
d  a 

2  hyperbolas 

(convex  branch) 

Figures  1,  2,  and  3  illustrate  the  elliptic  cases,  Figures  4  and  5  illustrate 
the  hyperbolic  cases,  and  Figure  6  illustrates  the  parabolic  case. 

Parabolic  and  hyperbolic  solutions  are  of  interest  for  a)  trajectories 
to  the  moon,  b)  some  satellite  intercept  problems,  and  c)  trajectories  of 
some  comets  in  the  solar  system. 

III.  DERIVATION  OF  LAMBERT'S  EQUATION 

Lambert's  equation  may  be  derived  directly  from  Kepler's  equation. 
The  basic  formulae  and  substitutions  involved  are  as  follows:^ 


r  =  a(l  -  e  cos  F)  (1) 

r  ,  +  r  =  2a  1  1  -  e  cos  4-  (F  +  F.)  cos  4  (F  -  F,)~!  (2) 

da  i-  2  a  d  2  a  d  J 

or,  if 


2G  =  F  +  F,  ] 
a  d 


and 


(3) 
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h* - 2o  - «h 

l — - - - - - — ( 

h*-2a-rQ  — 
U - 2 a-rfl  - H 


Figure  1.  Elliptic  Case:  4a  >  +  c. 


9 


i 


Page  9 


then 

Similarly, 

2 

c 

from 

Now  let 

so  that 


and  let 


then 


+  ra  =  2a  (1  -  e  cos  G  cos  g)  , 

4a2  sin2  G  sin2  g  +  4a2  (1  -  e2)  cos2  G  sin2  g 


2 

c 


2 

rd 


2r  ,  r  cos  0  . 

d  a 


cos  h  =  e  cos  G 


c  =  2a  sin  g  sin  h 
rcj  *■  ra  =  2a(l  .  cos  g  cos  h) 


c  =  h  +  g  ] 


> 


cos 


2  <« 


+  5)  = 


e  cos 


2  (Fa 


Fd> 


*  ra  +  c  =  4a  sin2  c 

r  .  +  r  -  c  -  4a  sin2  4  6 
da  Z 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


(10) 


(11) 


(12) 

(13) 
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and  substituting  in  Kepler's  equation 

nt,  -  F  -  F.  -  e(sin  F  -  sin  F,)  (14) 

f  a  d  a  d 

we  obtain 

ntj-  =  («  -  sin  t)  -  (6  -  sin  5)  (15) 

which  is  Lambert's  equation. 

Since  we  may  assume  that  F  -  F,  <2ir,  it  follows  that  0  <  (l/2)(c  -  6) 

a  d 

v  tt  and  0  (l/2)(i  +  6)<w,  so  that  0  <  (1 /2)  (  <  tr  and  -(1/2)  tt 

<  (1/2)S  <  (l/2)ir.  Hence  Lambert's  equation  has  four  possible  solutions; 
that  is,  combinations  of  each  of  the  two  solutions  for  t  and  6.  Let  t  j 
and  be  their  smallest  positive  values.  Then  the  four  solutions  are: 


Case  1A 

ntf 

=  <‘l 

-  sin  t  j)  -  (6  - 

sin  6  j ) 

(16) 

Case  2 A 

ntf 

=  (‘l 

-  sin  i  , )  t  (6j- 

sin  6  j ) 

(17) 

Case  IB 

ntf 

=  2rr 

-  (‘  j  -  sin  t  j)  - 

(6j  -  sin  6j) 

(18) 

Case  2B 

ntf 

II 

-  (t  j  -  sin  c  j)  t 

(6j  -  sin  6  j ) 

(19) 

which  correspond  to  the  four  possible  times  of  transfer  between  the  two 
points  in  the  ellipses  of  Figure  1. 

Call  0  the  angle  measured  counterclockwise  from  r^  to  rg.  Then 
Table  2  indicates  whether  the  motion  is  direct  or  retrograde. 


Table  2.  Direction  of  Motion,  Elliptic. 


and  /if 


0  <  0  <  IT 


2w  >  0  >  tt 


Case  1A  Direct 

Case  2A  Retrograde 

Case  IB  Direct 

Case  2B  Retrograde 


Retrograde 

Direct 

Retrograde 

Direct 
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The  equations  corresponding  to  motion  along  the  concave  branch  of  a 
hyperbola  are  developed  in  the  same  fashion  except  for  the  use  of  hyperbolic 
in  place  of  trigonometric  functions: 

r  =  a(e  cosh  F  -  1)  (20) 


(.  -  6)  =  <Fa  -  Fd) 


cosh  2  (*  f  M  =  e  cosh  (F^  +  F^) 


and  finally 


r  ,  4-  r  +  c  =  4a  sinh  c 
da  Z 


r  ,  t  r  -  c  -  4a  sinh  -r  6 
d  a  Z 


ntr  =  e(sinh  F  -  sinh  F.)  -  (F  -  F,) 
f  a  d  a  d 


nt^  =  (sinh  i  -  t  )  -  (sinh  6-6) 


The  above  equations  show  that  *  is  always  positive.  Furthermore  if  the 
angle  0.  as  defined  above,  is  less  than  ir  then  6>  0  and  if  0  >w,  then 
6  <0.  Therefore  two  cases  exist  for  hyperbolic  motion: 


Case  1C  ntj.  =  (sinh  c  -  c )  -  (sinh  6-6) 
Case  2C  nt^.  =  (sinh  t  -  c )  t-  (sinh  6-6) 

and  Table  3  tabulates  the  direction  of  motion: 


Table  3.  Direction  of  Motion:  Hyperbolic. 


and /if  0  <  9  <  n 


2w  >0  >  w 


Case  1C  Direct 


Retrograde 


Case  2C  Retrograde  Direct 
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Since  motion  along  the  convex  branch  of  a  hyperbola  is  of  no  practical 
interest  it  may  properly  be  omitted. 

The  equation  for  motion  along  a  parabola  was  found  by  Euler,  and  for 
completeness  it  will  be  included  here.  Euler's  equation  follows  from  Baker's 
equation^  (the  counterpart  of  Kepler's  equation  for  a  parabola): 


'vfTr (t  -  t)  v  ,  1  „  3  v 

^jTT/T'  tan  2  f  T  tan  I 


(28) 


Also,  for  a  parabola, 


2  v 

r  =  q  sec  j  =  q 


(i  +  tan2  j) 


(29) 


so  that 


r  ,  +  r  =  q 
d  a  M 


(  2  Vd  .  2  Va\ 

(2  +  tan  —  +  tan  —j—  I 


(30) 


and  the  equation  for  the  chord  becomes 


2  ,  .  *2  „  2  f  Va  ‘  vd  ^ 

c  =  (rd  ♦  ra)  -  4  rd  ra  cos  ^  ^ - j 


(31) 


or 


(vd  -  va) 


2^rd  ra"  cos  - =  H(rd  +  ra  +  c,,rd  +  ra  '  c)  <32) 


v  ,  v  a  (r  .  +  r  +  c)(r  ,  +  r  -  c) 

,  d  a  ,  V  d  a  d  a 

1  +  tan  -j"  tan  “J"  r  ^  1  -  1  """  ’  """  "" 


(33) 


so  that 


(rd  +  ra  +  c)  +  (rd  +  rfl  -  c)  +  2<y(rd  +  rg  +  c)(rd  +  r&  -  c) 


■( 


2q 


v  v  ,\2 

a  ,  d\ 
tan  ~2~  -  tan  -j—  I 


(34) 
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or 


Trd  +  ra  +  c)  {rd  +  ra  '  c) 


-  tan  —  “  tan  -j— 


(35) 


Using  equation  28  to  find  the  time  of  flight  between  two  points  in  the  orbit 
gives 


V 


f  .  va  . 

TJ2  =  tan  x  -  ta 


q 

which  can  also  be  written  as, 

^Ttf 


vd  ,  If. 

n  T  *  T  [tan 


3  va  3 

T  '  tan 


'ft  q 


I7T  =  ( 


tan  t*  -  tan 


-£■  j  3  ^1  +  tan  -p  tan  +  ^tan 


a  d 

~T  "  tan 


v 

c 

T, 


Substitution  from  equations  33  and  35  yields 

.3  /2  —  .  A  ,3/2 


6^tf  =  (rd  +  ra  +  c)  '  +  (rd  +  ra  -  c) 


(36) 


2n 


(37) 


(38) 


where  for  direct  orbits  the  upper  sign.  Case  ID,  is  used  if  0  <"  it  and  the 
bottom  one,  Case  2D,  if  0  >  rr.  The  signs  are  reversed  for  retrograde  orbits. 


Table  4.  Direction  of  Motion:  Parabolic. 

and  /if  6  <  tt  8  >  it 

Case  ID  Direct  Retrograde 

Case  2D  Retrograde  Direct 

IV.  APPLICATIONS 

As  previously  mentioned,  primary  applications  of  Lambert's  equations 
are  in  preliminary  orbit  determination  and  parametric  studies  of  transfer 
orbits.  This  application  involves  the  solution  of  transcendental  equations 
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and  the  problem  is  further  complicated  by  the  many  equations  which  are 

involved  (equations  16  through  19.  2 6,  27,  and  38).  The  most  practical 

4 

method  for  choice  of  equations  to  be  solved  is  presented  by  Breakwell 
and  involves  a  test  on  time. 


Breakwell's  parameters  are:  a  unitless  time 


T 


2  it  t 


f 


P 

s 


(39) 


where  Pg  =  period  of  elliptic  orbit  with  semimajor  axis  s/2;  a  unitless 
energy 


E 


energy  of  transfer  orbit 
E 

s 


(40) 


where  Eg  =  energy  of  elliptic  orbit  with  semimajor  axis  s/2;  and  a  unitless 
linear  scale 

K  =  1  -  -  (41) 

s 

where 

s  =  J  (rd  +  ra  +  c)  (42) 


A  plot  of  E  versus  T  for  a  K*>0.8  is- schematized  in  Figure  7.  When 
E  =  0  and  E  =  -1,  expressions  for  T  in  terms  of  K  are  easily  obtained 
and  become  landmarks: 


1A 

=  |  (1  -  k3/z) 

(43) 

2A 

=  |  (1  +  k3/z) 

(44) 

IB 

=  it  -  2  sin‘1^K,+  2 -  K' 

(45) 

2B 

=  it  +  2  sin"  I'^K'  -  2^1T^1  -  K1 

(46) 
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A  study  of  Figure  7  leads  to  Table  5  which  together  with  Tables  1  through 
4  allows  the  choice  of  the  proper  equation  to  be  used  and  determines,  with¬ 
out  further  solving  the  problem,  the  direction  of  motion  for  each  case. 


Table  5.  Classification  of  Cases  According  to  Time. 


If  tj.  is 


T2A  T1B 


T2A<T1B 


<T1A 

Then  Use  Case 

1C 

2C 

=  tia 

ID  (parabola  -  no  "solution"  nece 

>T2B 

IB 

2  B 

T2A 

2D  (parabola  -  no  "solution"  nece 

'  >T!A- 

<tib 

1 A 

2C 

*  T1B’ 

<T2A 

IB 

2C 

'  T2A- 

<T2B 

2A 

IB 

r>TlA' 

<T2A 

1A 

2C 

«  >  T2 A  ’ 

<tib 

1A 

2  A 

‘  >  tib’ 

<T2B 

2  A 

IB 


The  case  of  nearly  parabolic  motion,  either  from  the  elliptical  or  hyperbolic 
side  presents  some  difficulties.  Such  cases  have  been  solved  by  Lancaster^ 
where  the  elliptic  (or  hyperbolic)  equations  are  replaced  by  series  expansions 
which  converge  rapidly  for  nearly  parabolic  motion: 
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Case  1NP  T  = 


-  1 


(-1).  1.  3  •  5 •  •  •  (2n  -  1)  (-E)r 


n  =0 


->n-2  , 

2  n. 


2n  +  3 


n+3  /2 

(1  -  K  ) 


(47) 


Case  2NP  T  = 


i 

n  =0 


(-!)•  1  •  3 •  5 •  •  •  (2n  -  1)  (-E) 


n 


,n-2  , 

2  nl 


2n  +  3 


n+3/2 
(1  +  K  ) 


(48) 


Table  4  will  s-till  give  the  direction  of  motion,  where  Case  1NP  would  replace 
Case  ID  and  Case  2NP  would  replace  Case  2D. 

Equations  (47)  and  (48)  could  be  inverted  but  they  can  be  used  with 
iteration  as  they  stand  to  solve  for  E  when  T  is  given.  The  previous  tables 
are  summarized  in  Tables  (>  and  ~ . 


Table  h.  Summary  of  Cases. 


If _  Then  Use  Case  Equation 


‘f*T!A*  *T1 

1C  (26) 

2C  (27) 

T1A*  ATl<lf"TlA+  *T2 

1NP  (47) 

T2A-  AT3Stf  <-T2A+  *T4 

2NP  (48) 

T2A>T1B 

T1A  +  AT2  <lfS  T1B 

1A  (16) 

2C  (26) 

TIB<«f<T2A-  ATJ 

IB  (18) 

2C  (27) 

WT4<lf'T2B 

(17) 

IB  (18) 

T2A<T1B 

TlA  +  AT2<tf<T2A-^T3 

1A  (16) 

2C  (27) 

T .  ♦  AT  <  t  T 

*A  ai4  1  IB 

1A  (16) 

2A  (17) 

T  <  t  <  T 

IB  lf  1 2B 

2A  (17) 

IB  (18) 

lf  >  T2B 

IB  (18) 

2B  (19) 
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Table  7. 

Summary  of  Direction 

of  Motion. 

and  /if 

Cp  <  1 80 

cr  >  180 

Case  1C 

Direct 

Ret  rog  rade 

2C 

Ret  rograde 

Direct 

1  NP 

Di  rect 

Retrograde 

2NP 

Retrograde 

Direct 

1  A 

Di  rect 

Ret  rograde 

2  A 

Retrograde 

Di  rect 

1  AB 

Di  rect 

Ret  rograde 

2  AB 

Ret  rograde 

Di rect 

1  B 

Di  rect 

Ret  rog  rade 

2B 

Retrograde 

Di  rect 

Expressions  for  AT  j ,  AT^.  AT^,  ;ind  AT^  are  not  yet  available.  In 
practice  an  iterative  scheme  can  be  used  which  would  "jump"  from  Cases 
NP  to  C  or  NP  to  A  and  proceed  to  the  answer. 

Another  case  which  may  not  be  of  practical  use  is  when  E  approaches 
zero  from  the  IB  or  ZB  expressions;  that  is,  T  approaches  infinity.  These 
are  taken  care  of  by  the  following  series: 


Case  IBP 


Case  2BP 


T 


T 


2*  +  '■  (-!)•  1-  3-5-  .  •  (2n  -  1)  (-E)n 

(aR)3  L  l*"1  n!  2"  "  3 

2tt  \  (-  1  )•  1  •  3-  5.  .  .  (2n  -  1)  (-E)n 

(^fl)3  1  Z^ni  2n  +  3 


(1  +  Kn+3/2) 

(49) 


(1  -  Kn+3/2) 

(50) 


A  relatively  simple  computer  program  can  be  developed  which  would 
solve  all  the  above  cases  in  a  unified  iterative  procedure. 

The  above  study  is  limited  to  one  revolution  of  the  transfer  orbit.  For 
more  than  one  revolution  all  that  is  required  is  to  add  2Ntr  to  the  right  hand 
side  of  the  equations  to  be  solved  (where  N  =  number  of  complete  revolutions 
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of  the  transfer  orbit).  This  may  bring  in  some  additional  numerical 
difficulties  in  the  solutions  which  in  turn  may  produce  more  special  cases. 
Such  cases  have  yet  to  be  studied. 

V.  ORBITAL  ELEMENTS 


Once  "a"  has  been  obtained  from  the  solution  of  Lambert's  equations, 
the  rest  of  the  orbital  elements  may  be  easily  found/  For  elliptic  orbits 


ea  sin  F 


d 


(rfl  -  a)  -  (rd  -  a)  cos  (c  -  6) 
sin  (t  -  6) 


(51) 


ea  cos 


a 


which  gives  the  eccentricity.  Then 


M 


=  Fd- 


e  sin  F , 


where 


(52) 


(53) 

(54) 

(55) 


thus  obtaining  the  time  of  perifocal  passage.  To  find  the  angular  elements, 
the  cartesian  coordinates  or  the  direction  cosines  of  the  two  positions  must 
be  known. 


cos  i 


L2  +  L2  +  L2 

x  y  7. 


L 

x 


+  L2  +  L2 

y  z 


(56) 


sin  i  sin  0 


(57) 
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-  sin  i  cos  P  = 


L1  f  L2  f  l} 
x  y  v. 


(58) 


where 


-  y<  z?  -  y\  y 


1  '  2 
1  7'2 

V.  =  Xi  ^2  -  ^1  X2 


1  2 
1  x2 
;i  ^2 


(59) 


and  0  s  i  s  ir /2  for  direct  orbits  and  -nr / 2  c  i  tt  for  retrograde  orbits. 
Finally,  the  argument  of  perifocus,  w,  is: 


"  =  U1  *  v, 


(60) 


—  cos  v.  -  cos  F,  -  c 
a  d  d 


(61) 


~  sin  vd  =  \  I  1  -  e2  |  sin  F£ 


(62) 


where  and  may  be  obtained  from 


r^  cos  cos  0  +  sin  P 


r^  sin  u^  =  -x^  cos  i  sin  P+  cos  i  cos  P  +  z^  sin  i 


(63) 


For  hyperbolic  orbits,  the  following  substitutions  will  suffice 

1.  Use  hyperbolic  functions  for  the  eccentric  anomaly  (remember 
that  sin  -  sinh) 

2.  M.  =  -F.  t  e  sinh  F, 

d  d  d 

For  parabolic  orbits 


a  -  oo 
e  =  1 
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t 


and  for  time  of  perifocal  passage  use  Baker's  equation: 


. .  ,  Vd  ,  1  .  3  vd 

M  ,  =  tan  — r-  +  t  tan  -r- 

d  2  3  2 


where 


tan 


Av 

rd 

COS  -y-  ±  , 

r 

\ 

a 

sin 


Av 


Av 

-COS  — V-  ± 


tan 


T 


s  T  *  ^ 


sin 


Av 


(o  ^  ^ 


(Av  =  +0  for  direct  orbits 
Av  =  -0  for  retrograde  orbits) 


and  choose  v  ,  such  that  v  -  v  ,  =  Av.  Thus,  obtain  q  from 
dad 

rd  =  q  (l  ♦  tan2 

and  finally,  t,  from 


(64) 


(65) 


(66) 


The  angular  elements  are  obtained  as  for  ellipses  and  hyperbolas. 


(67) 


* 
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nomenclature: 

a  seminiajor  axis 

c  magnitude  of  cluird  between  two  radii 

e  eeeenl  rieit  y 

E  energy  ratio 

F  eeeent  r  ic  a  noma  1  y 

i  ini'lination  of  orbit  plane  to  referenee  plane 

K  see  equal  ion  4  1 

M  mean  anomaly 

n  mean  motion 

P  period 

q  radius  to  peri  foe  us 

r  magnitude  of  radius  vector 

S  semiperimeter 

t  time 

T  time  ratio 

v  true  anomaly 

h  variable  in  Lambert's  Equation  (see  equations  1  1  and  11) 

<  variable  in  Lambert's  Equation  (see  equations  li  and  li) 

fl  central  angle  between  the  two  radii  always  measured  counterclockwise 

from  r  ,  to  r  regardless  of  direction  of  motion 
d  a 

p  gravitation  and  mass  fac  tor 

t  time  of  perifocal  passage 

e  longitude  (or  right  ascension)  of  the  node 

w  argument  of  peri  focus 

Subindex 

d  departure 

a  or  rival 

f  flight 
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